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Abstract 

This paper is a tutorial that demonstrates various methods from the 
Colombeau theory of generalized functions in the context of semilinear 
wave equations. The Colombeau generalized functions constitute differ- 
ential algebras that contain the space of distributions. We solve the 1D- 
semilinear wave equation in these algebras, show how delta waves can be 
computed and then turn to linear and nonlinear regularity theory. 
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1 Introduction 

This tutorial serves the purpose of demonstrating how the theory of Colombeau 
algebras can be used to solve and study nonlinear partial differential equations 
that do not have solutions in the sense of distributions. It starts out with 
existence results obtained in the eighties of the twentieth century and leads up 
to recent regularity results. 

For the purpose of exposition we focus on a model problem, the one-dimensional 
semilinear wave equation with singular initial data 

d?u{x,t) -d 2 x u{x,t) = f(u(x,t)) + h(x,t), (x,t) eM 2 , , . 

u(x, 0) = a(x), d t u{x, 0) = b(x), xeR, ' 

where / is a smooth, polynomially bounded function and a, b and h are Colom- 
beau generalized functions on the real line and on the (x, t)-plane, respectively. 
The main questions to be addressed are the following: 

(a) existence and uniqueness of generalized solutions in Colombeau algebras; 

(b) limiting behavior of the representatives when the data are distributions; 

(c) regularity of generalized solutions. 

Concerning question (a), a short survey of existence and uniqueness results can 
be found in [25]. The existence of Colombeau solutions to semilinear hyper- 
bolic systems was one of the early applications of the theory [3J [23] . For lin- 
ear hyperbolic systems with non-smooth (Colombeau generalized) coefficients, 
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existence and uniqueness was established in }24| . for symmetric hyperbolic sys- 
tems in higher space dimensions in [18] and for hyperbolic pseudodifferential 
systems with Colombeau symbols in 15! . Notably for semilinear hyperbolic 
systems, question (b) has been answered in many cases, involving the notion of 
delta- waves QH E2 G6j E2 EH El] . Finally, regularity theory for Colombeau 
solutions is now based on the subalgebra Q°° of regular Colombeau functions 
and currently an active area of research, making use of pseudodifferential and 
microlocal techniques [pj [TT] [TTJ [2D] ■ In particular, the propagation of the Q°°- 
wave front set in linear systems with Colombeau coefficients is a theme of recent 
investigations [12]. To date, only few results are available about regularity of 
Colombeau solutions in the nonlinear case. 

The plan of the tutorial is as follows. After recalling the required notions from 
Colombeau theory in Section E] we shall prove existence and uniqueness of a 
solution u belonging to the Colombeau algebra on R 2 in Section El and state 
various preliminary regularity properties. Section [?] is devoted to computing 
the associated distribution (the distributional limit of the representing nets) 
when the initial data are delta functions or derivatives of delta functions. In 
Section E] we turn to regularity theory. We recall the <7°°-regularity result for 
the linear case and show that it fails in the nonlinear case. The tutorial finishes 
with a recent result on propagation of regularity in the nonlinear case in Section 

El 

For simplicity of presentation, we restrict our attention to the one-dimensional 
case and Lipschitz-continuous nonlinearity /. At the appropriate places of the 
tutorial, we will indicate what is known about the non-Lipschitz and the higher 
dimensional case. With Lipschitz-continuous /, in particular, the solution to 
([T]) exists globally in space and time. Therefore, we need not enter the discus- 
sion of domains of existence nor use other versions of Colombeau algebras to 
accommodate energy estimates when / grows polynomially. 

2 Notation 

The paper is placed in the framework of algebras of generalized functions intro- 
duced by Colombeau in [SJ [7] . We shall fix the notation and introduce a number 
of known as well as new classes of generalized functions here. For more details, 
see Q3]. 

Let fi be an open subset of M™. The basic objects of the theory as we use it are 
families (u e ) ee (o.i] of smooth functions u e 6 C°°(f2) for < e < 1. We single 
out the following subalgebras: 

Moderate families, denoted by Syi(Cl), are defined by the property: 

VK € £lVa e NS 3p > : sup \d a u £ (x)\ = as e -> 0. (2) 

xEK 

Null families, denoted by Af(Q), are defined by the property: 

VK € n\/a e N£Vg > : sup \d a u £ (x)\ = 0(e q ) as e -> 0. (3) 

xeK 

Thus moderate families satisfy a locally uniform polynomial estimate as e — > 0, 
together with all derivatives, while null functionals vanish faster than any power 
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of e in the same situation. The null families form a differential ideal in the 
collection of moderate families. The Colombeau algebra is the factor algebra 

= £ M (fl)/Af(n). 

The algebra Q(Q) just defined coincides with the special Colombeau algebra in 
[Ml Def. 1.2.2], where the notation G s (£l) has been employed. It was called the 
simplified Colombeau algebra in [4]. 

Families (r e ) e g(o,ii of complex numbers such that \r £ \ — 0(e~ p ) as e — > for 
some p > are called moderate, those for which |r e | = 0(e q ) for every g > are 
termed negligible. The ring C of Colombeau generalized numbers is obtained 
by factoring moderate families of^complex numbers with respect to negligible 
families. When SI is connected, C coincides with the ring of constants in the 
differential algebra G(fl). 

The restriction of an element u G Q{M. 2 ) to the line {t = 0} is defined on 
representatives by 

u\{ t =o} = class of (u e (-,0)) eE (o,i]. 

Similarly, restrictions of the elements of Q(fi) to open subsets of £1 are defined 
on representatives. One can see that SI — > 9 (CI) is a sheaf of differential algebras 
on 1". The space of compactly supported distributions is imbedded in Q(Ci) by 
convolution: 

l : £'(Q) -> <?(fi), = class of (to * (Ve)|n) e e(o,i]> ( 4 ) 

where 

<p e (x)=e- n <p(x/e) (5) 

is obtained by scaling a fixed test function ip G 5(]R n ) of integral one with all 
moments vanishing. By the sheaf property, this can be extended in a unique 
way to an imbedding of the space of distributions T>'(d). In the case SI = M n , 
the imbedding is given by the explicit formula 

i : P'(S1) -» Q(R n ), l(w) = class of (w * ( X <p B ))ee(o,i], (6) 

where x 1S some compactly supported smooth function identically equal to one in 
a neighborhood of zero. We refer to [IH Sec. 1.2] for further explicit expressions 
in general and in special cases. 

One of the main features of the Colombeau construction is the fact that this 
imbedding renders C°°(S1) a faithful subalgebra. In fact, given / G C°°(S1), 
one can define a corresponding element of C/(S1) by the constant imbedding 
er(/) = class of [(e,x) — > f(x)]. Then the important equality t(f) = <r(f) 
holds in 9(Q). 

If u G (?(S1) and / is a smooth function which is of at most polynomial growth at 
infinity, together with all its derivatives, the superposition f(u) is a well-defined 
element of Q(Q). 

The algebra G(p) can be equipped with a topology that turns it into a topo- 
logical ring and a topological module over C, and into a complete ultrametric 
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space. Let K g ft, m G No and p > 0. Then V(K,m,p) C £/(fi) is defined as 
the collection of elements of u G with a representative u e such that 

sup sup \d a u £ (x)\ — 0(e~ p ) as e — » 0. 

kG-R" |a|<m 

The sets V(.K', c 5(S1) define a base of neighborhoods of zero for a topology 
with the properties mentioned above. This topology was introduced in [3]; 
its vast potential was discovered by [39] who coined the term sharp topology. 
Important further developments are due to [TU] . 

We need a couple of further notions from the theory of Colombeau general- 
ized functions. Regularity theory is based on the subalgebra £7°°(f2) of regular 
generalized functions in Q(Q). It is defined by those elements which have a 
representative satisfying 

VK ^{l3p>0\/aeN% : sup \d a u e (x)\ = 0{e- p ) as e -» 0. (7) 

Observe the change of quantifiers with respect to formula Q ; locally, all deriva- 
tives of a regular generalized function have the same order of growth in e > 0. 
One has that (see Thm. 5.2]) 

^(fijnp'p) = c°°(n). 

For the purpose of describing the regularity of Colombeau generalized functions, 
G°°(Q) plays the same role as C°°(r£) does in the setting of distributions. How- 
ever, £?°°(f2) is not invariant under nonlinear maps. For this reason, various 
subspaces of G(Q) have been introduced to measuring regularity in the nonlin- 
ear case (see e.g. [29] ). We shall make use of just one of them, based on the 
notion of subsheaf regularity introduced in |19) . 

Definition 1. Cg(fl) is the space of the elements of Q(Q) with a representative 
(tt e ) e6 (o,i] such that lim e ^o u £ exists in L 1 1 oc (ri). 

Proposition 2. /(jCg(fi)J C Cg(fl) for every smooth function f all whose 
derivatives grow at most polynomially at infinity and which is Lipschitz contin- 
uous with a global Lipschitz constant. 

Proof: The polynomial bounds guarantee that f(u) is a well-defined element of 
Q(Q). Assume that u e converges to an element w G ^(K) on some compact 
set K. The estimate \f(u e ) — f{w)\ < Lipj|u e — w\ shows that f(u e ) converges 
to f(w), as desired. □ 

The kind of regularity that is encapsulated in the subspace above is described 
by the following rather obvious assertion: 

4(n)nx>'(n) = Li >c (n). 

We end this section by recalling the association relation on the Colombeau alge- 
bra Q(fl). It identifies elements of <?(f2) if they coincide in the weak limit. That 
is, u, v G G(Q) are called associated, u f» v, if lim £ _»o / \Ue(x) — v e (x))ip(x) dx — 
for all test functions if) G T>(Q). We shall also say that u is associated with a 
distribution w if u s —> w in the sense of distributions as e — > 0. 
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3 Existence/uniqueness of generalized solutions 

This section is devoted to solving the semilinear wave equation fl} in the 
Colombeau algebra Q(M. 2 ). Recall first that if w is a classical solution of the 
linear wave equation 

d 2 w(x,t) -d£w(x,t) = h(x,t), (i,i)el ! , . . 

w(x,0) = a(x), d t w(x,Q) = b(x), iel, { > 

then it solves the integral equation 
w{x,t) = -(a(x-t)+a(x+t))+~ b(y)dy+- / h(y,s)dyds. (9) 

Z Z Jx-t 1 JO J x — t+s 

Let Ko = [-K, k] be a compact interval. For < T < k, the trapezoidal region 
Kt is defined by 

K T = {(x, i) e R 2 : < * < T, \x\ < k - t}. (10) 
Using the following estimate is easily deduced (0 < t < T < n): 

||uj||l«(k t ) < l|a|| L ~(if ) + T \\ b \\h^(K ) + T ||/i|| L »o (jRrs) ds. (11) 

Jo 

For later reference, we note that estimate (fTTj) holds with the L 1 -norm in place 
of the L°°-norm as well. We now turn to the semilinear wave equation (fT]). We 
assume throughout that u — * f(u) is a smooth function all whose derivatives 
are of at most polynomial growth as \u\ — > oo, that / satisfies a global Lipschitz 
estimate (i.e., has a bounded first derivative) and that /(0) = 0. 

Proposition 3. Assume that the function f is as described above. Let a, b £ 
Q(R) and h G C?(R 2 ). Then problem {!]) has a unique solution u G Q(M. 2 ). The 
solution depends continuously on the data with respect to the sharp topology. 

Proof: To prove the existence of a solution, take representatives a e ,b £ ,h £ of 
a, b, h, respectively, and let u e G C°°(R 2 ) be the unique classical solution to the 
semilinear wave equation with regularized data: 

d 2 u £ - d 2 x u E = f(u e ) + h e on R 2 , 

u B (-,0) = o e , d t u e (-,0) = b e on R. 

The classical solution u e to (fT^j) is constructed by rewriting (TT2"|) as an integral 
equation and invoking a fixed point argument (this involves applying estimate 
(fTTj) successively to all derivatives). If we show that the net (u e ) ee ( ,i] belongs 
to £m(R 2 ), its equivalence class in £(R 2 ) will be a solution. We shall do the 
proof only for the upper half-plane (x, t) G R x [0, oo); the arguments for t < 
are similar. To show that the zero-th derivative of u £ satisfies the estimate ©, 
we take a region Kt and invoke inequality (fTT|) to see that 

IKI|L°°(K t ) < \\ae\\L°°(K ) +T\\b £ \\ L ^ {Ko) + T [ \\f{u e ) + h £ || L -(i<r s ) ds. (13) 

Jo 
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The last term on the right hand side of (fT5|) is estimated by 




Using that each of the terms involving a e ,b e ,h e is of order 0(e~ p ) for some 
p, we infer from Gronwall's inequality that the same is true of ||u e ||L°«(.R- t ) for 
< t < T. Thus u £ is moderate on the region Kt, that is, it satisfies the 
estimate ([2]) there. To get the estimates for the higher order derivatives, one 
just differentiates the equation and employs the same arguments inductively, 
using that the lower order terms are already known to be moderate from the 
previous steps. 

To prove uniqueness, we consider representatives u £ ,v £ £ £m[R 2 ) of two solu- 
tions u and v. Their difference satisfies 

df{u £ - v £ ) - dl(u e - v e ) = (f(u e ) - f(v e )) + n e , 
(u e (-, 0) - v £ (-, 0)) = n 0e , d t (u £ (-, 0) - v £ (-, 0)) = n le 

for certain null elements n £ ,nQ £ ,ni £ . Thus u £ — v e satisfies an estimate of the 
form (|13p . but with the null elements n £ , uq £ , n\ £ replacing a £ ,b £ , h £ there. This 
implies as above that the L°°-norm of u £ — v £ on Kt is of order 0(e q ) for every 
q > 0. By [TU Thm. 1.2.3], the null estimate ([3|) on u £ — v £ suffices to have null 
estimates on all derivatives. Thus u = v in (?(IR 2 ). 

The proof of continuous dependence on the data a,b,h follows exactly the same 
lines. □ 

Having established a general existence- and uniqueness result, the question arises 
what else can be said about the solution, when more is known about the data. 
A number of results on the qualitative properties are indeed available. We begin 
with the most basic result on C °° -smoothness of the solution when the data are 
smooth. 

Proposition 4. Assume that a,b belong to C°°(IR) and h is in C°°(IR 2 ). Then 
the generalized solution u £ ^(R 2 ) to problem (QP coincides with the classical 
solution w G C°°(R 2 ) ; that is, u = l{w) in G(R 2 ). 

Proof: The imbeddings l and a coincide on C°°(R) and on C°°(M 2 ). Thus we 
may represent the data by a £ = a, b £ = b and h £ = h. The proof of Prop. [3] 
shows that u £ = w is a representative of the generalized solution, and this means 
that u = l(w) in G(R 2 ). □ 

Proposition 5. Assume that the data a, b and h are continuous functions and 
let w G C(R 2 ) be the corresponding continuous (weak) solution. Then the gen- 
eralized solution u G £/(IR 2 ) is associated with w, that is, u « l(w). 

Proof: Let a £ — a * tp £ and similarly for b £ and h £ with appropriate one- or 
two-dimensional mollifiers (p £ . On the one hand, the families a £ ,b £ ,h £ define 
representatives of the generalized functions t(a), i(h). On the other hand, 
they converge locally uniformly to the original continuous functions a, b, h. The 
estimates (fT"3|) and (flT]). applied to u £ — w, a £ — a etc. together with Gronwall's 
inequality, show that the solution u £ converge to w with respect the sup-norm 
on compact sets. In particular, u is associated with w. □ 
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Note that in general u ^ l(w) as elements of t/(M™ +1 ). As a rule, equality of 
continuous and generalized solutions does not hold in any algebra of generalized 
functions. This is similar to Schwartz' impossibility result and has been worked 
out in [25 . When the data are distributions, there may be no meaning for 
a distributional solution, in general. Yet the solution in Q(K. 2 ) may still be 
associated with a distribution. We will turn to an incident of such a situation 
in the next section. 

Proposition [3] is a model result. In fact, for vanishing driving term h = 0, it is a 
special case of [3] and of [23] : for Lipschitz continuous, smooth /, existence and 
uniqueness of a solution in Q(R 2 ) can be proven in space dimensions n = 1,2,3, 
see e.g. [28]. If / is not Lipschitz, but of polynomial growth, energy estimates 
can be used to construct solutions in the Colombeau algebra 52,2(8" x [0, 00)) 
introduced in [5]. As in the classical case, the growth type of / is connected 
with the space dimension; the cases 1 < n < 6 have been treated in [21] : for 
n = 3 see also [TJ. 



4 Delta waves 

In this section we investigate the behavior of the solution u € Q(M. 2 ) to the 
semilinear wave equation (|TJ) when the initial data a, b are given by distributions. 
In certain cases, the smooth solutions u s with regularized initial data a e ,b £ 
converge to a distributional limit. When it exists, this limit is called a delta 
wave. Stated equivalently, the generalized solution u 6 (?(R 2 ) with data t(a), t(b) 
admits an associated distribution. To be sure, this distribution is not a solution 
of the original nonlinear wave equation in any sense. It rather has to be viewed 
as describing the qualitative properties of the Colombeau solution. 

For simplicity of exposition, we will assume here that the right hand side h 
is zero (see e.g. [28] for a result with nonzero h) and will exhibit the model 
case of bounded nonlinearity / and initial data with discrete support. Thus 
a, b G 2?'(R) will be assumed to be Dirac measures or their derivatives situated 
at a finite set of points. Consider the weak solution v S 2?'(R 2 ) to the linear 
problem 



d 2 v - d 2 x v = on 
v(-,0) = a, dtv(-,0) — b on 



(15) 



By d'Alembert's formula, it is given by 



1 / , s , ,x 1 



x+t 



v(x,t) = -(a(x-t)+a(x + t)) + -J b(y)dy. (16) 

Its singular support S consists of the characteristic lines with slope ±1 emanat- 
ing from the points of support of a, 6, hence has Lebesgue measure zero. Denote 
by I5CW the almost everywhere defined, measurable function equal to v(x,t) 
for (x, t) g" S. In fact, the function l s c v is piecewise constant and arises from 
the measure part in b solely. For example, when b(x) = S(x), we have that 
l s cv(x,t) = 1/2 for \x\ < t and lscv(x,t) = otherwise. If the function / is 
continuous, then the nonlinear equation 

dtw - d 2 w = f(l s c v + w) on R 2 , 

w(-,0) = 0, d t w(-,0) = on K { ' 
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has a unique weak solution w G C(R 2 ). Finally, let u G £(R 2 ) be the Colombeau 
solution to the nonlinear problem 

d 2 u-d 2 u = f(u) on R 2 , 

u(-,0) = i(a), d t u(-,0) = i(b) on M. 1 ' 

As we shall see shortly, v + w is a delta wave, i.e. the associated distribution 
corresponding to the Colombeau solution to (1181) . 

Proposition 6. Assume that the data a,b G 2?'(R) have discrete support, f is 
smooth, bounded and globally Lipschitz. Then the generalized solution u G (?(R 2 ) 
is associated with the sum of the distributions v G 2?'(R 2 ) and the continuous 
function w G C(R 2 ), that is, u w l(v + w). 

Proof: The solution u G 5(R 2 ) has a representative which satisfies 

d 2 u £ - d 2 u £ — f(u e ) on R 2 , 

u £ (-,0) = a e , d t u £ (-,0) = b £ on R 

where a e = a * (p £ ,b £ = b * ip e . Further, let v £ be the classical smooth solution 
to the linear wave equation (|15[) with data a £ ,b £ . We have 



(dfu - dl){u £ -v e -w) = f(u £ ) - f(l s c v + w) on 
(u £ — v £ — w)(-, 0) = 0, d t (u £ — v £ — w)(-, 0) = on 



(19) 



We rewrite the right hand side of the first line in (|19|l as 

f(u £ ) - f(v £ + w) + f(v £ +w)- fOsc v + w). 
Using (1 1 1 p with the iZ-norms, we arrive at 



\\u e ~ V E - to||Li(Kt) < T \\f'\\L°°(R)\\Uc ~ V £ - w\\ L i( Ks )ds 

Jo 

+ T 2 \\f(v £ +w)- /(l s c v + w)\\ lHKt) 

valid for < t < T. Since f(v £ + w) — f(ls c v + w) converges to zero almost 
everywhere and remains bounded, Lebesgue's theorem shows that its L^-norm 
on Kt converges to zero. By Gronwall's lemma, it follows that the L 1 -norm of 
u £ — v £ — w converges to zero on any Kt as well. Hence u £ converges to v + w 
weakly, which translates into the claimed association result. □ 

This result is an example of a wide range of much more general results on 
the existence of delta waves. We refer to the papers [331 123 1221 13Z]- In case 
the driving term h is white noise, the generalized solution to the semilinear 
(stochastic) wave equation is associated with the solution of a linear stochastic 
wave equation in many cases. This has been shown e.g. in [TJ [301 White 
noise or positive noise in the initial data has been considered in |31[ 135] . 



5 Linear regularity theory 

Consider the linear wave equation (| 1 5|) with initial data distributions a, b G 
2?'(R). The classical propagation of singularities results says that the singular 
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support of the solution v £ 2?'(R 2 ) is contained in the union of characteristic 
lines emanating from the singular support of the initial data. For example, when 
the singular support of the initial data is situated at the origin, the solution is a 
smooth function outside the forward and backward light-cone S = {(x, t) £ R 2 : 
|x| = \t\}. In the Colombeau setting, this is not so, as shown by the following 
example. 

Example 7. Let u £ £?(R 2 ) be the Colombeau solution to the linear wave equa- 
tion 

d?u-d 2 u = Q on R 2 , 

u (.,0) = 0, d t u{;0) = l(S) 2 on R. 

with initial data given by the square of the Dirac measure in Q (R) . Then u does 
not coincide with a C°° -function in the interior of the light-cone {(x,t) £ R 2 : 
|x| < \t\}. 

In fact, by d'Alembert's formula, the solution u £ C/(R 2 ) has the value 

1 f°° (\ f°° 

— J t((5) 2 (x)dx — class of (- J Lp 2 (x)dx 

inside the forward light-cone {(x, t) £ R 2 : |x| < t}, t > 0, which is a generalized 
constant not belonging to C. 

The example shows that the space C°° is not suitable for measuring regularity 
in Colombeau algebras. Its role is taken over by the algebra G°° defined in the 
Introduction. We have the following result. 

Proposition 8. Let v £ £/(R 2 ) be the Colombeau solution to the linear wave 
equation H5\) with initial data a, b, £ (?(R). If a,b £ (J°°(R \ {0}) ; then u £ 

g°°(m\s). 

Proof: We take a cut-off function \ S C°°(R) which is equal to one in a small 
neighborhood of the origin and zero otherwise. Then v = v± + vi where v\ £ 
t/(R 2 ) is the solution with initial data a\ = (1— %)a, b\ = (1— x)b and V2 £ £?(R 2 ) 
is the solution with initial data 02 = x a i ^2 = xb- Using d'Alembert's formula, 
we have that 

dxVu = ^(a l£ (x + t) + a' le (x - t) + b le (x + t) - b u (x - t)) 

and similarly for dtV\ s . Thus if all the derivatives of a e and b e are of the same 
local order, so are the derivatives of v\ e . This shows that v-y belongs to Q°°(M. 2 ). 
Similarly, differentiating v^ once shows that the first and higher derivatives of 
V2 vanish outside a small neighborhood of S, determined by the diameter of the 
support of x- Since x was arbitrary, this shows that V2 is a generalized constant 
on R \ S, thus belongs to g°°(R \ S) as well. □ 

It is clear that Prop.[S] can be generalized to initial data with discrete Q°°- 
singular support. It also holds in any space dimension 27J. 
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6 Nonlinear regularity theory 

This section addresses propagation of singularities for the semilinear wave equa- 
tion. The initial data that are assumed to be regular outside the origin, but 
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may have a singularity there. Following the classical theory of propagation of 
jump discontinuities one may hope to prove that, in the one-dimensional case, 
the solution will be regular (in the appropriate sense) in the interior of the 
light-cone S. That this behavior occurs in the C°°-category for the semilinear 
wave equation in one space dimension was shown by [38]; it does not hold in 
higher space dimensions or for higher order operators |36j . for which anoma- 
lous singularities may occur. The question we ask is whether this transport of 
regularity into the interior of the light-cone happens in the Colombeau setting 
as well. We do not have a general result yet, but we will be able to prove this 
transport of regularity under conditions which resemble the situation leading to 
delta waves. We shall employ the subalgebra G°° and the space Cg defined in 
the Introduction for measuring the regularity. 

We consider the semilinear wave equation |T]). As before, we assume that the 
function u — > /(it) is smooth with all derivatives of at most polynomial growth 
as \u\ — > oo, and that it satisfies a global Lipschitz estimate, i.e., has a globally 
bounded derivative. We take initial data of the form 

Ui = n + Si e g(R), i = o,i 

where 

r, G £g(R) and supp(s 4 ) = {0}, i = 0, 1. 

Define the generalized complex number MeCas the class of 

1 f°° 
M £ = - / s u (x)dx 

" J — oo 

where (sie) E e(o,i] 1S a representative of S\. We shall write M w m if M e con- 
verges to a complex number m G C and \M\ « oo if \M £ \ — > oo as e — * 0. Again 
S = {(x,t) € K 2 : \x\ — \t\} denotes the forward and backward light-cone. 

Proposition 9. Assume that the function f is smooth, globally Lipschitz, and 
all its derivatives are polynomially bounded. Let u £ ^(R 2 ) be the solution to 
the semilinear wave equation flp with initial data uq,u\ as described above. If 
either 

(a) f is globally bounded and M « m for some mGC, or 

(b) f is globally bounded, hin^^oo f(y) exists and \M\ ~ oo 
then u 6 g°°(M. 2 \S) +C 1 g (R 2 ). 

Proof: We first derive the estimates for (x, t) in the upper half-plane R x [0, oo). 
Denote by Is the characteristic function of the solid light-cone S = {(x, t) : t > 
0, |x| < t}. Let (rj e ) £e (o 1 i], (sie) e e(o,i]) * = 0, 1, be representatives of the initial 
data, where we may assume that supp(si £ ) c [—1], rj\ with rj as small as we wish. 
Let u e , v £ e C°°(M x [0, oo)), w £ G L 1 1 oc (R x [0, oo)) be the solutions to 

(9f - d%) u e = f(u e ), u e (-,0)=u Oe , d t u e (-, 0) = u le , 

(9 2 - dl) v e = 0, w e (.,0) = SOe, d t v s (;0) = s u , 

(d 2 t -dl)w e = f{M e lv+w £ ), W£ (.,0) = r Q£ , d t w e {-,0) = r le . 

Then 

(9 2 - d 2 x ) (u £ -V £ - We ) = f(u e ) - f(v £ + W £ ) + f(v £ + W £ ) - /(M E 1 E + W E ) 
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with zero initial data. The difference of the first two terms on the right hand 
side can be estimated by the L°°-norm of /' times \u £ — v s — w e \, while the 
difference of the last two terms vanishes off an 77-neighborhood of the light-cone 
S, because v e (x,t) = M £ for \x\ < t ~ r\. Taking a trapezoidal region Kt as 
in (fT0|) . the boundedness of /, inequality (jTTJ) and Gronwall's lemma give an 
estimate of the form 

IK - v £ - w e \\ h i {KT) < Cn (20) 

for some constant C > and all e 6 (0, 1]. In the case (a), let w € L 1 1 oc (R x 
[0, 00)) be the solution to 

{pi - dl) w = /(ml s + w), w(-,0)=r o , d t w{-,0)=r 1 

where Vi is the limit in L 1 1 oc (R) of ri £ as e — > 0, i = 0, 1. Now 

(dl -dl) (w e -w) = f(M e l s +w e )-f{M e lv + w) + f{M e l s +w)-f(mlv+w) 

with initial data given by ri £ — i = 0,1. The difference of the first two 
terms on the right hand side is bounded by the L°°-norm of /' times \w £ — w\. 
By assumption, the difference of the last two terms converges to zero almost 
everywhere. Using inequality (jTTJ) with iZ-norms in place of the L°°-norms and 
Gronwall's lemma as above shows that 

\\w e - HIlH-Kt) ase^O. (21) 

In case (b), defining w S Ll oc (M. x [0, 00)) as the solution to 

(Of - dl) w = Ll s , tu(.,0)=r , d t w(;0)=n, 

the same argument as above leads to the convergence result (|21[) in this case 
as well. Combining ([2TJ]) which holds for arbitrarily chosen r\ > with (|2"T| 
shows that u £ — v £ converges to w in 1}(Kt) as e — > 0. By Prop. [51 v £ enjoys 
the ^°°-estimate (O off the light-cone S. Thus u £ = v £ + (u £ — v £ ) defines an 
element of 0°°((R x [0,oo))\5) +C 1 g (Rx [0,oo)). 

The arguments for (x, t) in the lower half-plane M x (— 00, 0] are the same with 
—M in place of M. □ 

Remark 10. The hypotheses on the generalized constant M in are 
satisfied when the term s± in the initial data is a polynomial in the Dirac measure 
and its derivatives. In fact, when s± — d a i(6), we have M w 1 for a = and 
M w for a > 0. When s\ = n(i(5)) for some polynomial function n, only the 
cases M « m for some to G C or \M\ sa 00 can occur. 

Prop. shows that regularity of the type Q°° + £g is propagated into the region 
inside the light-cone. It is clear that it can be generalized in various ways: for 
example, the support of the singular part sq,si of the data could consist of a 
discrete set rather than just a point. However, it should be noted that continuous 
dependence of the regularized solutions on the data in terms the L^-norm enters 
into the argument, and such a property depends decisively on the particular 
equation, the space dimension and the nonlinearity /. Prop. exploits such 
special properties and thus falls short of providing a prototypical description of 
nonlinear propagation of regularity for Colombeau solutions, which remains a 
challenging open question. 
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Another promising approach to measuring regularity in the nonlinear case is 
the Colombeau-Holder-Zygmund-scale which has been introduced and applied 
to nonlinear scalar first order equations in [16] . 
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